Disorder induced transverse delocalisation in ropes of carbon nanotubes. 
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A rope of carbon nanotubes is constituted of an array of parallel single wall nanotubes with nearly 
identical diameters. In most cases the individual nanotubes within a rope have different helicities 
and 1/3 of them are metallic. In the absence of disorder within the tubes, the intertube electronic 
transfer is negligeable because of the longitudinal wave vector mismatch between neighboring tubes 
of different helicities. The rope can then be considered as a number of parallel independent ballis- 
tic nanotubes. On the other hand, the presence of disorder within the tubes favors the intertube 
electronic transfer. This is first shown using a very simple model where disorder is treated perturba- 
tively inspired by the work in reference-. We then present numerical simulations on a tight binding 
model of a rope. Disorder induced transverse delocalisation shows up as a spectacular increase 
of the sensitivity to the transverse boundary conditions in the presence of small disorder. This is 
accompanied by an increase of the longitudinal localisation length. Implications on the nature of 
electronic transport within a rope of carbon nanotubes are discussed. 

PACS numbers: 



I. INTRODUCTION 

A rope of single wall carbon nanotubes (SWNT) is 
generally made of ordered parallel tubes with different 
helicities, but with a narrow distribution of diameters^. 
The center of the tubes form a triangular lattice so that 
there is for each metallic tube in a rope on average two 
neighboring tubes which are also metallic. In the ab- 
sence of disorder within the tubes, the intertube elec- 
tronic transfer, defined as the matrix element of the 
transverse coupling between two neighboring tubes, in- 
tegrated over spatial coordinates, is negligeable because 
of the longitudinal wave vector mismatch between tubes 
of different helicities^. The rope can then be consid- 
ered as made of parallel independent nanotubes. The 
transport is ballistic and the dimension less conductance 
(in units of the quantum conductance Gq = 2e 2 /h) is 
two times the number of metallic tubes within the rope. 
However, it has been shown^ that disorder within the 
tubes favors intertube scattering by relaxing the strict or- 
thogonality between the longitudinal components of the 
wave functions. Using a very simple model where disor- 
der is treated perturbatively, we show in section II that 
the intertube scattering time is shorter than the elas- 
tic scattering time within a single tube. In tubes longer 
than the elastic mean free path, this intertube scatter- 
ing can provide additional conducting paths to electrons 
which would otherwise be localized in isolated tubes. 
In the limit of localised transport along the tubes we 
show that the longitudinal localisation length is not a 
monotonous function of disorder and increases at mod- 
erate disorder. In order to go beyond these analytical 
results we have performed numerical simulations on a 
tight binding model of coupled ID chains with different 
longitudinal hopping energies. This model described in 
section III mimics the physics of transport in a rope of 
carbon nanotubes in the sense that in the absence of dis- 
order the electronic motion is localised within each chain. 
Transverse delocalisation as a function of disorder is in- 



vestigated through the sensitivity of eigen-energies to a 
change of transverse boundary conditions from periodic 
to antiperiodic. 

These results show that disordered ropes of carbon 
nanotubes can be considered as anisotropic diffusive con- 
ductors, which in contrast to individual tubes, exhibit a 
localization length that can be much greater than the 
elastic mean free path. 



II. ELECTRONIC STRUCTURE OF ROPES OF 
CARBON NANOTUBES 

A. Band structure of a rope without disorder 

We consider a rope constituted from SWNT with di- 
ameters ranging between 1.2 and 1.5 nm. It can be shown 
(see table HJ that the tubes within such a rope can have 
different kind of helicities. Following the model devel- 
oped by Maarouf and Kane^ one can characterize the 
electron wave functions at the Fermi energy e p with two 
wave vectors k± and k\\ respectively perpendicular and 
parallel to the tube axis : \^ >= \k± > \k\\ > such that 
ty(x,y) ~ e lkj - v e lk w x . From this model, it is possible to 
compute the matrix elements of the transverse coupling 
hamiltonian H±_ between 2 tubes a and b: 



(9 a \H ± \Hf b )=t±{a,b)5 kaltkH (1) 

with t ± (a,b) = f Te -(V4)flao(fca-L-fci,-0 2 aiK f t T = 
tG\/-£jfR — 7.5 meV for an average tube radius R = 
0.7 nm, tc — 0.1 eV is the inter-plane hopping energy 
in graphite and ao = 0.5 x 10~ 10 m is the Bohr radius. 
The term Sk all ,fc b n is due to the othogonality between wave 
functions of different longitudinal wave vectors and shows 
that tubes of different helicities are uncoupled to first or- 
der in H±. The exponential part takes into account the 
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FIG. 1: Top, Brillouin zones of two metallic tubes of dif- 
ferent helicities. Bottom: Dispersion relation as function of 
the longitudinal component of the wave vector of the low- 
est energy bands for these 2 tubes for which 7^ kZi 
at the Fermi energy illustrating the second order coupling 
|-ff±| = 2t±(a, b) 2 /8e(a, b) via the higher energy state 
at energy cf + 5e(a,b)/2 for which k a \\ = fcjji. 



mismatch between the tt orbitals belonging to neighbor- 
ing tubes of different helicities 3 . 

Going beyond first order perturbation it is possible to 
show that neighbouring tubes of different helicities are 
coupled to second order in H±_ with a characteristic en- 
ergy scale ej_2 = 2t±(a,b) 2 /5e (Fig. [T]) involving transi- 
tions to higher energy states for which e b (k) = e a (k) = 
(F + 8e/a,b)/2. This second order coupling gives rise 
to an intertube hopping probability which is very small 
compared to the inverse ballistic time L/vf of a micron 
long nanotube. It is thus reasonable to assume that, in a 
rope constituted of tubes of different helicities, electronic 
transport is confined within each ballistic tube and ex- 
hibits a strong ID character. 



B. Disordered ropes in the perturbative regime 

In the presence of disorder, plane waves localised on 
individual tubes are perturbed into wave packets and 
become much more sensitive to the transverse coupling 
leading to a transport regime which is no longer ID 
but can be delocalised on percolating clusters of metallic 
tubes within the rope. 

We consider a very short range (on site) disorder poten- 
tial, such that its average < V(x) >— and its variance 



y/< V(x) 2 > = W. 



(2) 



where the index a runs on the chains constituting the 
rope, each of them being characterized by its atomic sites 
x a and the the disorder potential W(x a ) which probabil- 
ity distribution is given by: 



P(W(x a )) = 1/W if w(x a ) £ [-W/2- W/2] 
otherwise 

. (3) 

The disorder perturbed wave functions can be written 
to first order in disorder: 



I*' >=|fcio> [\kp>+Yl 



■\\o) 



e(0) -e(i) 



\k\u > 



(4) 



where i runs over unoccupied states and V(k) is the 
Fourier component of V at wave vector k. We do not 
consider transitions involving different values of k±_ and 
neglect possible perturbation of t±(a, b) with disorder. 

This perturbed wave function contains plane waves of 
all values of parallel momentum. As a result electrons 
can hop from tube to tube and conserve their momentum. 
The intertube coupling energy between 2 tubes a and b 
is now E±_ 

the disorder potential 



< v f , nl#jJ v f , b > which reads to first order in 



2 x-^ W(x a )cos((k\\ a — k\\h)x a ) 
E ± (a,b) =< k a± \H ± \k b± >-J2 (h " 77 T 

(5) 

The disorder average value of this coupling is zero but its 
typical value e± = y E± is equal to: 



e±(a,b) = t ± (a,b) 



W 



3<5e(a, b) 



(6) 



with 6c(a, b) = Hvf \k\\ a — fc|if>|- 

This disorder induced intertube coupling energy is re- 
lated to the second order coupling term calculated in 
the previous section through: e±(a, b, W)/e±2{a, b) — 
W/t±(a,b) which, as will be shown below, can be much 
larger than one in a typical rope of SWNT. The coef- 
ficients tx(a,b) and |^7^y were calculated for all the 
helicities corresponding to metallic tubes given in table U 
with diameter between 1.2 and 1.5 nm, and are depicted 
in the histogram fig. ([2]). From these values we obtain 
the average value e± = 0.03 W . 

To investigate the nature of transport in ropes, it is in- 
teresting to compare the typical intertube hopping time 
T~h = ft/ (e_i_) to the intra tube scattering time r e in- 
duced by the same disorder. The related elastic mean 
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helicity (n,m) 


(9,9) 


(10,10) 


(11,8) 


(12,9) 


(12,6) 


(13,7) 


(15,6) 


diameter (rim) 


1.24 


1.38 


1.32 


1.45 


1.26 


1.4 


1.5 


helicity (n,m) 


(14,5) 


(13,4) 


(16,4) 


(15,3) 


(17,2) 


(16,1) 


(18,0) 


diameter (nm) 


1.36 


1.23 


1.46 


1.33 


1.44 


1.32 


1.43 



TABLE I: Values of possible helicities and diameters of 
the 14 metallic tubes with a diameter between 1.5 and 1.2 
nm. Diameter and helicity are related through: D(n, m) = 
2^ V m 2 + ran + n 2 . There are also 24 insulating tubes in 
the same diameter range. 
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FIG. 2: Histogram of possible values of the coefficients t j_ (a, b) 
and for all couple of helicities (a, b) of metallic tubes 

with a diameter between 1.2 and 1.5nm. The average yields 
3e bot = 0.09. 



free path l e — vpT e was calculated by White et Todorov 9 

and found to be given by l e — j^knc for a tube of nc car- 
bon atoms along the circumference, where tp = 2.7 eV 
is the Fermi energy (measured from the bottom of the 
band) and W 2 is the variance of the disorder (assumed 
to be short range). This value of l e is unusually large 
compared to what is expected in an ordinary conductor, 
since it is proportional to the number of sites along the 
circumference of the tube. This is due to the existence of 
only 2 conduction modes at the Fermi energy regardless 
of the diameter of the tube. As a result there exists a 
rather large range of disorder W for which r e is greater 
than Th — h/ (ex), which means that a charge carrier 
can visit several neighboring tubes between 2 elastic col- 
lisions. As shown schematically on fig[3] comparing the 
relevant time scales for transport in a rope of micron 
length, four different regimes (1,2,3,4) can be reached as 
the amplitude of disorder is increased. 

1. At very low amplitude of disorder, when both r e 
and Th are long compared to the ballistic time t = 
L/vf, transport is ballistic and one dimensional and 
electronic wave functions are localised within single 
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FIG. 3: Different transport regimes depending on the ampli- 
tude of disorder compared in a rope of length L = lfi m as 
discussed in the text. Region 1 = ID ballistic, region 2 = 3D 
ballistic and region 3 = 3D diffusive. 



tubes. 

2. When Th < r& < r e transport is still ballistic but 
wave functions are delocalised over several tubes. 

3. When r& < Th < T e transport is diffusive along 
the percolating clusters of metallic tubes along the 
rope as long as the rope is shorter than the local- 
isation length which typical value is £ — 2N m l e ; 
where N m is the number of metallic tubes in the 
largest percolating cluster of metallic tubes within 
the rope. Typical values of N m for commonly in- 
vestigated SWNT ropes in experiments are given 
in the appendix. 

4. At large disorder m, < r e < Th electronic states are 
localised within individual tubes at the scale of l e 
(not shown in fig[3]). 

From this qualitative model one expects that trans- 
verse transport but possibly also longitudinal transport is 
favored when increasing disorder in the rope. We present 
in the following section numerical simulations which con- 
firm this statement. 



C. Analytical results in the localised regime 

In the following we consider the case where, in the 
absence of inter chain coupling, electronic wave functions 
are localised within each tube aligned along the x axis 
and can be characterised by the set of parameters x a , k a 
and the localisation length £ a such that: 

^ a (x) = cos(k a (x - x a ) exp-(|.T - X a \/£ a )/y/(£ = 

^h^Q L (fa(k ~ k a ) + fa{k + K) exp [£((x - X a )] 

. (?) 

where f a is a Lorentzian function centered on k a of width 
5k a = 27r/c; a . In the presence of a small intertube cou- 
pling such as described by eq. 1121 one can easily compute 
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the typical transverse coupling energy between two near- 
est neighbor tubes at lowest order in t± : 



<* a |flx|*6 



>= 



exp ik(x a — X(,)X 4 



(8) 



XL 



(\ 2 + (k±k a ) 2 )(X 2 + (k±k b ) 2 ) 



In the following we assume that 8k a — Skb = A = 27r/£ 
the summation runs on multiple values of 2tt/L up 
to n = L/ao the number of sites on the tubes of length 
L. Averaging over the random phase factors expi(k — 
k')(x a — xt,) for k ^ k' leads to the average square of this 
coupling: 



T ab =\<^ a \H ± \^ b >\ 2 

4^(a,6) 2 /(AL) 2 E fe 



A 



[A 2 + ((fc±A :a ) 2 ] 2 [A 2 + ((fc±fc b ) 2 ] 2 

(9) 

Keeping in the summation over k only the \/(2tt/L) 
terms centered around k a and kb within 8k = A, finally 
yields: 



T ab = 4(t ± (a,b)) 2 \L 



A 2 



ttL 2 (A 2 + (k a - fc,) 2 ) 2 



(10) 



A, the inverse typical localisation length for ID disor- 
dered chains can be approximated by A = Ca^W 2 / E F 
where the amplitude W of the intra chain on site disor- 
der is assumed to be small compared to the Fermi energy 
Ep. One can see from expression (JTUJ) that the typical 
transverse coupling between tubes of different helicities 
k a =h kb obtained after averaging on the positions x a and 
Xb increases with disorder like W 6 at low disorder and de- 
creases like 1 /W 2 at large disorder such that A > \k a —kb\. 



coupling becomes zero to first order in t± since identi- 
cal longitudinal Bloch wave vectors correspond to differ- 
ent energies i ra cosfcao, except at half filling kao = tt/2. 
Note however that this disorder among the t n also im- 
plies a distribution of Fermi velocities which does not 
exist in CNT. This leads to the following hamiltonian 



H = H 



N m N s -1 
n=l n s = l 



H \\=J2Y1 tn n s ><n,n s + l\ + h.c] (11) 



*.lE£=i[e 



In, n s >< n + 1, n s | + 



exp(i(j))\N m , n s >< 1, n s \ + h.c. 



(12) 



The last term in H±_ corresponds to the periodic 
boundary conditions around the cylinder modeling the 
rope. This periodic boundary condition involves a phase 
factor exp(i<fi) equivalent to a fictitious flux $ through 
the cylinder which modifies the phase <fi = 27r$/$o of 
transverse boundary conditions of the wave functions 6 . 
Due to the different values of t n the wave functions are 
localised within each chain in the limit of very long chains 
N s 3> N m and t± <C tn. An extra disorder hamiltonian 
Hd is added either as a random distribution of on-site 
potentials Wi of width W or as an extra random con- 
tribution St n (s,s + 1) to the nearest neighbor coupling 
t n within the chain n (bond disorder) characterized by a 
distribution of width 6t w assumed to be independent of 
n. 



Numerical results. 



III. A SIMPLIFIED ANISOTROPIC TIGHT 
BINDING MODEL FOR A ROPE OF SWNT. 

We investigate in the following a simplified tight bind- 
ing model for the transport in a rope containing a perco- 
lating cluster of N m metallic SWNT. As shown in the ap- 
pendix, we expect that in commonly investigated SWNT 
ropes N m is at most equal to 10. Each SWNT labeled n 
is described by a simple ID atomic chain of N s sites with 
a nearest neighbor coupling energy along the chain t n 
which can be different from one chain to the other. The 
variables t n are randomly distributed around their aver- 
age value fii with a square distribution of width 8t\\. For 
convenience we take the chains on the surface of a cylin- 
der where only sites belonging to nearest neighbor chains 
are coupled (see FigH]). This crudely reproduces the sit- 
uation of a hexagonally packed rope, where each metallic 
tube has on average 2 metallic tubes as nearest neigh- 
bors. The interchain transverse coupling is described by 
a transverse nearest neighbor coupling t± <C tu. The dis- 
tribution among the t n plays the same role as the helicity 
distribution among the tubes in a rope i.e. the inter chain 
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FIG. 4: Left: Schematic modelisation of a rope of nanotubes 
by N m slightly coupled tight binding chains of different longi- 
tudinal energies ti 7^ t% 7^ £3 etc.... Right: The chains lye on 
a cylinder threaded by a fictitious Aharonov Bohm flux which 
enables to investigate transverse transport. 



In order to investigate the interchain localisation and 
the effect of disorder along the chains, we have calcu- 
lated the eigenvalues of H(<j)) = Hq(4>) + Hd- As shown 
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in the seminal work of Thoulessi the sensitivity of these 
eigenvalues to a change of the phase <f> of the boundary 
conditions can be considered as a measure of the delo- 
calisation of the wave functions on the various chains 
constituting the rope. More precisely we have computed 
the quantity Se perp —< \ek(4> = i") — efe(0)| >k where 
the average <>k is taken on the N m N s /A energy levels 
of the spectrum between 1/8 and 7/8 filling excluding 
the region between 3/8 and 5/8 filling. Around half fill- 
ing the tight binding dispersion relations for all chains 
are indeed crossing each other whatever the values of t n 
are and the model is inadequate. For ropes of N m = 10 
chains of N s = 100 sites the phase dependent eigenvalues 
were determined using a standard matlab diagonalisation 
routine. On the other hand for longer ropes (N s = 1000) 
a Lancsos algorithm^ was used to compute the eigenval- 
ues spectrum. The quantity Se perp calculated for ropes 
of ten tubes is shown on Figj5] as a function of disorder 
strength both for on site and bond disorder. 

When the t n are all equal (5tt\ = 0), Se perp is a 
monotonously decreasing function of disorder as expected 
in the physics of standard localisation. On the other 
hand when the t n are different, Se perp is very small at 
low values of disorder as expected, with a power law de- 
pendence in (this exponent corresponds to the hop- 
ping probability around a circumference with 10 sites), 
see FigjBJ More important, at fixed value of tj_, 5e perp 
increases with disorder amplitude, goes through a maxi- 
mum and decreases at large disorder. As shown in Fig [5] 
the conditions of observation of this non monotonous de- 
pendence of Se perp as a function of disorder amplitude 
depends drastically on the amplitude of the transverse 
hopping integral. It is clearly observed for very small 
values of t±_ with a low disorder increase in \y 2 - 5±0 - 5 fol- 
lowed by a decrease in 1/W 1 . This power law exponent 
is consistent with the analytical result derived in previous 

1 /2 

section for the quantity r Q ' b oc 1/W describing the large 
disorder hoping between two adjacent tubes, its exten- 
sion to hopping processes around a rope containing N m 
chains yielding a decrease in l/WQV m ). When increasing 
t_i_ multiple order hopping processes in f j_° dominate the 
transverse transport and Se perp becomes independent of 
disorder at low value of W ■ 

The sensitivity to a phase shift along the chain di- 
rection was also investigated from the computation 
Se pa r —< |efc( < / ) || = — efc (0) | >k where <j>\\ is the phase 
factor on periodic boundary conditions parallel to the 
tube axis. For long ropes (1000 sites along the longitu- 
dinal axis) it is possible to observe an increase of Se par 
with t±, see Fig. [7] This behavior is associated with 
an increase of the longitudinal localisation length with 
t± in the range of disorder where the inter tube cou- 
pling increases with disorder. One can easily deduce from 
this figure that for the value of W — 1.3 the localisation 
length increases by approximatively a factor 4 at 

t±/t\\ = 0.06. This is done assuming an exp(— L/£(t±)) 
behavior for 5e va r- 
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FIG. 5: Transverse delocalisation by disorder: Evolution of 
Stperp with the amplitude of disorder for 10 coupled chains 
of 100 sites, corresponding to W for on site disorder, and 
5t w for bond disorder. The situation of identical values of i n 
corresponding to St\\ — exhibits expected disorder induced 
localisation, whereas the situation with different values of t n 
where Stu/tu = 0.2 shows a regime of disorder induced trans- 
verse delocalisation . The transverse hopping energy is chosen 
to be tx = 0.06tii. v is the density of states (inverse nearest 
level spacing). 



IV. CONCLUSION: IMPLICATION FOR THE 
TRANSPORT IN ROPES ON CARBON 
NANOTUBES. 

We have shown that a rope of single wall carbon nan- 
otubes of different helicities is expected to exhibit very 
different regimes of electronic transport depending on 
the amount of disorder. When disorder is very small 
electronic transport takes place in independent ID bal- 
listic tubes. The conductance is then expected to be 
G = 2NGq where N is the number of connected metal- 
lic tubes. The experimental observation of a strong shot 
noise reduction in short low resistive ropes^ confirms the 
ballistic nature of transport in these ropes. On the other 
hand disordered tubes are expected to behave as 3D dif- 
fusive multi channel conductors whose maximum value 
of conductance is G = 2N m (l e / L)Gq when L is smaller 
than the localization length £ = 2N m l e , where N m is the 
number of metallic tubes in the largest percolating cluster 
of metallic tubes in the rope. Four probe transport mea- 
surements on NT ropes after ion irradiation damage^ 
have been shown to give information on the intertube 
hopping processes within the rope which was found to 
increase with disorder. These findings can also be inter- 
preted taking into account the increase of the intertube 
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FIG. 6: Results obtained by diagonalisation of a system 
of N m — 10 chains of N s = 1013 sites. Continuous 
curves, average energy difference 8e perp as a function of on 
site disorder amplitude. From top to bottom: tx/tii = 
0.1,0.07,0.03,0.02,0.01. Note the dependence in t^O at low 
disorder. Dashed curve power law fits approximation those 
curves with 2.5 and —9 exponents (left and right part). 



scattering rate with disorder. The physics of intertube 
transfer has also been shown to play an important role 
in the transport in multiwall nanotube o 11 ' 12 . The situa- 
tion is however different than in ropes since each tube is 
only coupled within first order to two other tubes (near- 
est inner and outer shells) and moreover only the most 
external tube is connected to electrodes. 

Let us also mention that a similar scenario of disorder 
induced delocalisation has been predicted to take place in 
networks of disordered polymers as discussed i n 13 ' 14 i 15 . 

We finally note that these different types of trans- 
port are also expected to influence the superconductiv- 
ity observed at very low temperature on ropes of carbon 
nanotubes. The superconductivity in weakly disordered 
ropes has been observed to exhibit a strongly ID char- 
acter with a T=0, H=0 transition. On the other hand 
more resistive ropes exhibit a broad transition at finite 
temperature characteristic of a multi channel quasi 3D 
syste m 16 ' 17 . 

Acknowledgments: We acknowledge very fruitful dis- 
cussions with Christophe Texier, Nicolas Dupuis, Michael 
Feigelmann, Piotr Chudzinski and Alexei Ioselevitch on 
this problem. 
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FIG. 7: Results obtained by diagonalisation of a sys- 
tem of N m = 10 chains of N s — 1013 sites. Aver- 
age typical energy difference 8e par between the symmet- 
ric/antisymmetric boundary conditions in the longitudinal di- 
rection as a function of t±/t\\. From top to bottom W/t\\ = 
0.55,0.9,1.1,1.3,1.5 



V. APPENDIX: DETERMINATION OF THE 
TYPICAL SIZE OF PERCOLATING CLUSTER 
OF METALLIC TUBES WITHIN A ROPE. 



Carbon nanotubes in a rope are arranged according 
to a triangular network, with on average 1/3 of metal- 
lic tubes and 2/3 of semiconducting ones. Since 1/3 is 
below the percolating threshold, 1/2, for nearest neigh- 
bor couplings in the 2D triangular lattice, the metallic 
tubes do not percolate over the full rope and constitute 
disconnected clusters which size distribution depends on 
the number of tubes within the rope. We have calculated 
numerically this size distribution for ropes containing a 
few hundred tubes. This result is shown on FigJS] for 
ropes containing 100 and 400 tubes. The size distribu- 
tion decays exponentially with the number of tubes in a 
cluster. We find that for a rope containing 400 tubes the 
average cluster size is 3 but by integrating the number of 
clusters of size above a given value we find that there is 
at least one cluster of tubes of size above 15. 
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FIG. 8: Size distribution of the number of tubes n in a per- 
colating cluster of metallic tubes for ropes containing N=100 
and N=400 tubes among which 1/3 on average are metallic. 
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